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Abstract 
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kind of conditions leads to the formation of naked singularity. For this purpose, we apply the 1+3 
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collapse of a perfect fluid. This formalism allows us to construct an autonomous system of evo- 
lution and constraint equations for scale-invariant dynamical variables normalized by the volume 
expansion rate of the timelike orthonormal frame vector. We investigate the asymptotic evolution 
of such dynamical variables towards the formation of a central singularity and present a conjecture 
that the steep spatial gradient for the normalized density function is a characteristic of the naked 
singularity formation. 
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I. INTRODUCTION 



The investigation of the end states of gravitational collapse of sufficiently massive stars is 
a long standing problem to be pursued in general relativity. One of the remarkable results is 
the naked singularity formation and this is one of the possible end states in spherically sym- 
metric gravitational collapse starting from regular initial data sets with nonzero measure. 
This was firstly shown by the detailed analysis of the Lemaitre-Toleman-Bondi (LTB) so- 
lution describing the inhomogeneous dust gravitational collapse, and has been subsequently 
confirmed in several other matter models including perfect fiuid models 
In relation to the cosmic censorship conjecture, much attention therefore has been paid to 



the generality and the physical reasonability of the initial data sets 



singularity formation in the various matter models (see 



Qy,y, 
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eading to the naked 



ll| for example). 



Besides the initial data analysis, dynamical aspects of spherically symmetric gravitational 
collapse at later stages have been studied in terms of the shear function associated with the 
velocity vectors of the collapsing matter near the center [l2, and it has been suggested 

that the shear can contribute to the delay of the apparent horizon formation. This role of 
the shear function may be important in considering the dynamical mechanism of the naked 
singularity formation. In this paper, however, we would like to focus on another dynamical 
approach to the end-state problem in spherically symmetric gravitational collapse. Our main 
concern is the relationship between asymptotic behavior of dynamical variables towards the 
formation of a central singularity and the causal structure of the arising singularity. 

As an useful method to analyze such dynamical properties, we adhere to the H-3 or- 
thonormal frame formalism which was originally proposed by Elst and Uggla et al. [l^ 
to study the dynamical behavior of the gravitational field variables such as the shear aab and 
the expansion B associated with the timelike orthonormal frame vectors near the spatially 
inhomogeneous cosmological initial singularity. This formalism is based on the coordinate 
independent representation of Einstein's field equations in the form of an autonomous system 
of the first order evolution equations and constraints with the scale-invariant dimensionless 
variables normalized by the Hubble scalar H = Q/3. By virtue of this Hubble normalization 
of dynamical variables, it is possible to remove the time dependent factors due to the volume 
contraction given by the rate 6 measured in the local reference frame. As will be shown 
in this paper, the Hubble normalized density function Q can approach zero or remain finite 
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with the lapse of time towards singularity formation in the spherically symmetric inhomo- 
geneous perfect fluid collapse, even though the central proper density becomes divergent at 
this final stage. 

Our main purpose in this paper is to discuss the relation between the asymptotic behavior 
of Q and the causal structure of the singularity. For the LTB solution, using the causal 
structure classified by the papers [l, 2], we find that the growth of the steep spatial gradient 
of the profile of fl near the center in the case — > is a characteristic property leading 
to the naked singularity formation. If the density function Q remains finite at the final 
stage, the end-state problem becomes more subtle. Nevertheless, we can discuss the critical 
value of the density contrast which gives a threshold of the transition from the black hole 
formation to the naked singularity formation. 

This paper is organized as follows: In Sec. [Tll we begin with a brief review of the 1+3 
orthonormal frame formalism with the Hubble normalized variables and apply it to a spher- 
ically symmetric perfect fluid system. In this formalism, we adopt the separable volume 
gauge, which specifies the lapse function to be equal to the inverse of the Hubble scalar 
H. In Sec. IIIIl considering the inhomogeneous dust gravitational collapse described by the 
marginally bound LTB solution, we show that this gauge condition is useful to examine the 
asymptotic behavior of the Hubble normalized variables towards the central singularity for- 
mation. Then, we relate the asymptotic behavior of the Hubble normalized density function 
Q to the arising causal structure of the end states of collapse. We propose a conjecture that 
the larger spatial gradient of the asymptotic profile of Q is essential to the naked singularity 
formation. In Sec. IIVI to support this conjecture, the asymptotic analysis is extended to 
gravitational collapse of perfect fluid with pressure. Taking account of the causal struc- 
ture of sp herically symmetric self-similar spacetimes which has been clarified by previous 
works jl?!. isl . [l9 |. we numerically estimate the critical value of the density contrast using 
the asymptotic profile of Q for the naked singularity formation. The results are summarized 
in Sec. |Vl Throughout this paper, the units in which SttG = c = 1 are used. 



II. BASIC EQUATIONS FOR HUBBLE NORMALIZED VARIABLES 



In this section, following the 1+3 orthonormal frame formalism developed in 



15 



we present an autonomous system of evolution equations and the constrains for the scale- 
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invariant variables in a spherically symmetric system with a perfect fluid source. We express 
an orthonormal frame as {cq, Sa} (where a = 1,2, 3) with the unit vectors Bq and 
representing the timelike reference congruence and the rest of 3-spaces, respectively. The 
frame metric is given by 7]^^, = diag[— 1, 1, 1, 1]. For simplicity, we specify the timelike frame 
vector Co to be hypersurface orthogonal, and the spacelike frame vectors to be nonrotating 
Fermi-propagated along the integral curves of Cq. 

Defining the four-velocity u = eo for the unit timelike vector Bq tangent to the reference 
congruence, we can introduce the basic geometrical quantities (-it", 6, through 
the commutator relations as follows 

[eo, Ba] = iiaeo - (^^f + , (2.1) 

[e„. Bp] = (2aiJpJ + e^/ssn^'^^ , (2.2) 

where the square brackets denote the antisymmetric part of a tensor, and e^/j^ is the totally 
antisymmetric three dimensional permutation tensor. The scalar function B is the volume 
expansion rate, the vector -u" and the tensor cr°^ (the trace-free symmetric tensor) are the 
acceleration rate and the shear rate of the frame vector eo, respectively. They are calculated 
from the equations 

e = vx , (2-3) 

iia = u'^Vf.Ua , (2.4) 

0-a(3 = ^ipUa) - ^ (r^a/S + UaUp) + U{aUp) , (2.5) 



where the round brackets denote the symmetric part of a tensor. (See [15|] for the definition 
of the covariant derivative in the orthonormal frame formalism.) In addition, the quan- 
tities a" and n"^ (the symmetric tensor) determine the connection on the one-parameter 
family spacelike hypersurfaces which can be defined by the assumption of the hypersurf ace- 
orthogonality of the timelike frame vector eo. 

Now let us turn our attention to the introduction of the basic variables characterizing the 
matter field. We consider a perfect fluid as the collapsing matter in this paper. We assume 
the equation of state to be 

p = {i- (2.6) 

with the constant 7 lying in the range 1 < 7 < 2. The pressure p and the energy density jl 
are measured by a comoving observer who has the same velocity as the fluid four-velocity u. 



In general, the fluid four-velocity u is not equal to the four-velocity u defined by the timelike 
frame vector cq, thus we introduce the basic matter variables fi and v by decomposing the 
energy-momentum tensor of the perfect fluid with respect to u into the form 

T^,u = lJ'{uf,u^ + lG~^ {2v(^^v^) + v<^v^y)} +p{g^^ + u^u^) , (2.7) 

where the angle brackets denote the symmetric trace-free part of a tensor. We have the 
following relations 

^ = r^G/i , p = |7 - 1 + (^1 - ^7) ^^X} (2-8) 

with the scalar functions G and F are deflned by 

G ^ 1 + (7 - l)vX , r = ^ . (2.9) 

The vector v represents the peculiar fluid velocity relative to the rest 3-spaces of eo, and 
defined through the relations 

with the fluid four- velocity u normalized as u^u'^ = —1. 



16| is to introduce the scale- 



An important procedure of the formalism developed in 
invariant dimensionless variables by normalizing the geometrical and matter variables using 
the Hubble scalar 

= (2.11) 
We denote the Hubble normalized quantities as 

^0 = 1 , 9^^^, (2.12) 
{f7",S,^,A-,iV"/5} = -i{M",a,^,a",n"^} , (2.13) 
^ = (2.14) 

where Q is the Hubble normalized density function. The vector v is a dimensionless vari- 
able and the Hubble normalization is not necessary for this variable. To introduce a local 
coordinate system, we also deflne the Hubble normalized components of the frame vectors 
as ^ 



As expressed in Eq. (12.151) . we hereafter attach the hat " to spacetime indices in order to 
distinguish them from the orthonormal frame indices. As physically interesting additional 
scale-invariant variables, we introduce the deceleration scalar q and the spatial Hubble gra- 
dient Aq, defined by 

q = -l-^doH, (2.16) 
\^ = -l-d^H. (2.17) 

12 

They will be used to eliminate the Hubble scalar H appearing in the evolution equations 
and the constraints.^ 

Some gauge choice is still allowed within the framework of the 1+3 orthonormal frame 
formalism, and in |l^, [l^, the evolution equations and the constraints for the Hubble nor- 
malized variables have been given with the so-called separable volume gauge, which simplifies 
the temporal frame derivative Oq to 

^0 = -dt (2.18) 
using a nondimensional time coordinate t."^ Further, the additional gauge constraint 

Ua = K (2.19) 



is required for the separable volume gauge (see for its details). Owing to this specification 
of the gauge, the matching of the four-velocity u = oi the reference congruence with 



the fluid four- velocity u (i.e., t> = 0) as in 20| is not always permitted. The important 



geometrical result of Eq. (12.181) with the commutator equation (12.11) is that the volume 
density V defined by V^^ = det(ej) has the form 



V = Vo X e-^\ (2.20) 

where Vq is an arbitrary function of spatial coordinates. From Eq. (I2.20p . in the limit t oo, 
the volume density approaches zero to form a singularity. This is a useful property of the 



^ Although the spatial Hubble gradient was expressed as in [16|, we have changed its notation in order 
to prevent readers from confusing it with the radial coordinate r which will be introduced later. 

^ While the lapse function is the positive definite function equal to the Hubble scalar H in the separable 
volume gauge of 16], in this paper we specify the lapse function to be —H in order to keep the positivity 
of the lapse function. The minus sign in the right hand side of Eq. p.lSp . which does not appear in 



the corresponding equation in 16|, comes from this gauge specification, which may be referred to as the 



separable volume gauge for gravitational collapse. 
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separable volume gauge to investigate the asymptotic dynamical behavior just before the 
singularity formation in gravitational collapse. 

Now let us study spherically symmetric gravitational collapse with perfect fluid using 
the Hubble normalized scale-invariant variables. By virtue of the Hubble normalization, 
the Hubble scalar H becomes the only variable carrying a physical reference scale and the 
analysis of the evolution equations and the constraints for these variables will allow us to 
observe dynamical behavior deviated from the time dependence of the volume contraction 
of the reference congruence. We consider the spherically symmetric line element of the form 

ds'^ = -I'^{t,r)dt^ + J'^{t,r)dr'^ + R^{t,r){de'^ + sm^ed(t?) . (2.21) 

The coordinate expressions for the orthonormal frame derivatives can be given by 

eo = r^di , Sr = J~^dr , ee = R'^dg , = sin'^ 98^ . (2.22) 

As is used in Eq. (12.221) . hereafter we substitute the letters r,6,(j) (instead of the numbers 
1,2,3) into the spacetime spatial indices i and the orthonormal frame spatial indices a (i.e., 
i = r,9,(l) and a = r,9,(f)). The gauge conditions (12.181) and (12.201) are reduced to 

/ = -H-^ (2.23) 

and 

J{t,r) = C{r)R-\-^\ (2.24) 

where C is an arbitrary function of r. 

An autonomous system of evolution equations and constraints for the Hubble normal- 
ized variables presented in 15|, ll6| has been expressed with the orthonormal frame deriva- 
tives applicable to any spacetimes. It is a straightforward to apply the formalism to 
the spherically symmetric metric (I2.2ip with the separable volume gauge. The quantities 
{E/ , Eg^ , , X'^ jT,"^^ ,Q,v^ , q) are independent Hubble normalized variables to be analyzed 
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here and we arrive at the following evolution equations for these variables 

E; = -{q- S'-'^) E/ , (2.25) 

e/ = - (^g + ^S-^ e/ , (2.26) 

A' = -{q- S"') A' - ^E/{E''y , (2.27) 

y = -{q- S") A'- - E/q' , (2.28) 

t'-' = -[q - 2)s^- + \e;{a')' - \E;{yy - ^x-a' - \{Es'f - 2-iG-^n{v'-f , (2.29) 

r / 9 \ 



n = -{2q - 1)^] + ?>G-' |7 - 1 + - -7J {vy^ n (2.30) 

+ 7E/(G~l^]^;'^)' + -fO'^nv'^ {v'J:'''' - 2A') , 
and the constraints 

1 + ^ ^2E/{Ay - 2AM'^ - 3{Ay + (V)'} " ^(S")' - = , (2.31) 
E/(S"'')' + 2A" - - 3A"S"' + 3-fG~^Qv' = , (2.32) 
E/{Efy - (A'- + y)E/ = , (2.33) 

where the dot and the prime mean the partial derivatives dt and dr, respectively. These 
equations correspond to the Einstein equations, the Jacobi identities and the contracted 
Bianchi identities. In addition to these equations, we can use the following Raychaudhuri 
equation for the deceleration scalar q: 



l + 3G-N7-H-(l-^7)(^T 



n . (2.34) 



It is remarkable that no time derivative of the radial velocity appears in these set of the 
evolution equations and the constraint equation (12.321) can be used to determine v^. We can 
also check that the remaining two constraints fl2.3ip and (12.331) are consistent with other six 
evolution equations for {E/, Eg^, A"^, X^, E^'^, Q). 

Finally, let us explicitly present the Hubble normalized variables using the metric func- 
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tions and the matter fields: 



E; = ~ , E,' = -L , (2.35) 
= -Z\l_ :i\ , (2.36) 



J ' 






' j 


-l[ 




IdrR 




JR 




dri 




J 





A' = , (2.37) 

y = , g = -1 - 7 , (2.38) 



u 



(2.39) 



3 {1 - (i-r) 

These relations are useful to analyze the asymptotic behavior of dynamical variables. Al- 
though the description of time evolution with the Hubble normalized variables and the 
separable volume gauge is available until the singularity is formed at t = oo, the spacetime 
region covered by the radial coordinate r may be too restricted to determine the nakedness of 
the arising singularity. Thus we will restrict our investigations to the models of gravitational 
collapse of which causal structure of the end states is already known. An example of such 
a model is the LTB solution, which is well known as a generic model of the inhomogeneous 
dust gravitational collapse. In the next section, we will examine the dynamical behavior of 
the Hubble normalized variables in the marginally bound LTB spacetime to discuss their 
key feature relevant to the naked singularity formation. 

III. ASYMPTOTIC BEHAVIOR IN INHOMOGENEOUS DUST COLLAPSE 
A. Brief review of the LTB solution 

Let us begin with a brief review of the LTB solution describing inhomogeneous dust 
gravitational collapse (see for its details). This solution includes the two arbitrary 

functions of radial coordinate usually denoted as F and / , which are related to the Misner- 
Sharp mass and the initial velocity, respectively. In this paper, we consider only the solution 
with / = 0, which is called the marginally bound solution. Using the comoving coordinate 
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system {r, p, 6, 0}, the line element for this solution is 



ds^ = -^dT^ + {dpRfdp^ + {d9^ + sin^ OdcP^) , (3.1) 

2/3 



R{r,p)=p{l~\l^{l + T)\ , (3.2) 



where B is an arbitrary positive constant. 

The time r = — 1 is usually interpreted as the initial time at which the area radius R 
becomes equal to the coordinate radius p. In addition, without loss of generality, we can 
choose r = to be the time of the central singularity formation by specifying the leading 
order term of the arbitrary function F{p) near the regular center p = as F ~ i?p^. In the 
comoving coordinate system, the components of the dust four-velocity are given by 

= , uP = u^ = ui' = 0, (3.3) 

and from the Einstein equations we obtain the proper energy density p measured by a 
comoving observer as 

The approximate form of the function F near the regular center p = can be written as 

F(p) = V (1 - 2F„p" + 0(p"+^)) (3.5) 

with a positive integer n. The subleading term F„p" in Eq (13.51) represents the dominant 
inhomogeneity of an initial dust distribution near the center, and we assume n > 2 and 
Fn > to require the regularity of the proper energy density at the center (i.e., dpp = 
at p = and at r = —1). The causal structure of the end states of the collapse is closely 
related to the inhomogeneity given by the term F^p^- For = 2, the shell- focusing naked 
singularity appears at the center p = at the time r = 0. For n > 4, the arising singularity 
is hidden behind the event horizon. For = 3, using the parameter 

(3-6) 



the condition for the naked singularity formation is given by 
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B. Asymptotic behavior of the Hubble normalized variables 



In this subsection, we would like to clarify how the asymptotic behavior of the Hubble 
normalized variables leading to the naked singularity formation depends on the initial den- 
sity inhomogeneity characterized by the integer n. With the help of the relation between 
the causal structure of the end states and the choice of n obtained in , we discuss what 
asymptotic behavior of the Hubble normalized variables characterize the causal structure of 
the singularity. As the formation of the central singularity occurs at the point r = p = 
in the comoving coordinate system {r, p, 6*, 0}, our strategy is to analyze the asymptotic 
i-dependence of the Hubble normalized variables in the limit t — oo by using the coordi- 
nate system {t, r, 6*, 0} with the separable volume gauge condition. For this purpose, we 
consider the coordinate transformation between the two coordinate systems, which leads to 
the following partial differential equations for r(t,r), p{t,r), I{t,r) and J{t,r) 

- {d,Rf(? = I' , (3.8) 
{d,Rfp" - = J' , (3.9) 

{d,Rfpp' = ^rT'. (3.10) 

We also have Eq. (12.241) as the separable volume gauge condition. As will be shown, the 
four Eqs. ( I3.8p - (]3.10p with ( 12.240 demand that timelike curves with r = const, to converge 
to the singular point r = p = 0. The behavior of these coordinates is schematically shown 
in Fig. [H 
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T 




FIG. 1: A schematic diagram describing the relation between the comoving coordinate systems 
{t, p} and the separable volume gauge coordinates {t,r}. In the limit t oo, all the timelike 
curves labeled r = const, (solid curves) converge to the onset point r = p = of the singularity 
formation. 

In particular, the exponential t- dependence of the function JR^ in Eq. f l2.24p significantly 
affects the asymptotic relation between the two coordinate systems {t, p} and {t,r}, and 
the converging behavior of the r = const, timelike curves is useful to analyze the asymptotic 
dynamical features just before the central singularity formation. In this subsection, we ana- 
lyze the asymptotic t-dependence of the Hubble normalized variables on such a congruence 
of the timelike curves. When the coordinate variables r(t, r) and p{t, r) approach zero in 
the limit t — > cxd with a fixed value of r, the metric function R{t, r) also goes to zero as 



The key issue to be analyzed here is which of the terms -F„p" and r in Eq. (13. lip becomes 
dominant in the limit t ^ oo along the timelike curves. 

From Eq. (I2.24p . we assume that the metric functions J and R have exponential t- 
dependence in the limit t ^ oo. This assumption turns out to be compatible with Eqs. (13. 8p - 
(I3.10p under the relations 



R{r, p) c:^ p{F^p" - t)'/' . 



(3.11) 



Rr^Jr^rr^Ir^ eXp(— t) 



(3.12) 
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if the exponential t- dependence of p is also derived from Eq. (13. lip . To check this, we 
first assume that the ratio F„p"/|r| approaches zero. In this case, Eq. (13.111) leads to 
p ~ exp(— i(:/3), which is consistent with the assumption F„p"/|r| only for n > 4. 
If the ratio -F„p'"/|r| is assumed to blow up, we have p ~ exp{— 3t/(3 + 2r2)}. It is easy 
to check that this is allowed only for n = 2, which gives the exponential t-dependence as 
p ~ exp(— 3t/7). For the remaining case n = 3, we obtain p ~ exp(— 1/3) and the ratio 
FnP^ /\t\ remains finite. 

Now let us see the asymptotic behavior of the Hubble normalized variables using the 
asymptotic exponential forms of the metric functions J, J, R and the coordinate variables 
r, p. Because we have the relation Eq. (13.121) irrespective of the choice of n, Eqs. ( ^1 



(I2.38P mean that the Hubble normalized variables {EJ, Eg^, A^, X^) become finite and their 
values can only depend on r. On the other hand, variables S'"'" and q approach zero. Using 
the condition -u^ = for the comoving coordinate system, the radial velocity given by 
(I2.39P can be rewritten into the form 

(-"f = ■ (3-13) 

Hence, the radial velocity remains finite in the limit t oo irrespective of n. From 
Eqs. (13:81)- (13101) ■ we have 

and using this relation, the Hubble normalized density Q given by Eq. ( 12.401) becomes 

4f2n2 

Q= S ■ (3.15) 
9R^dpR ^ ' 

Applying the approximate form ( 13. lip of R and the asymptotic relation f ~ — r, we have 
the following asymptotic form of the density function 

4 

^ 3{3 + (2n + 3)F„pVkl}{l + i^nPVkl} ' ^^'^^^ 
Owing to the term F„p"/|r| contained in Eq. (I3.16p . an interesting difference of the asymp- 
totic behavior of Vt appears according to the choice of n and this will be shown in the next 
subsection. 



13 



C. Relation between the asymptotic Hubble normalized density and the causal 
structure of the end state 

We rewrite the term F„p"/|r| in Eq. fl3.16p as a function of t and r. In the hmit t — > oo, 
the metric functions can be written as J = /o(r) exp(— t), J = JqIv) exp{—t) and R = 
Ro{r) exp{—t) . Because the choice of the spatial coordinate r remains arbitrary within 
the framework of the separable volume gauge, the arbitrary function C{r) is included in 
Eq. fl2.24p for J. To remove this ambiguity, we introduce the new spatial coordinate ( by 

-E/^ = l, (3.17) 

where the Hubble normalized variable E/ = —Iq/Jq should be regarded as a function of r. 
This specification of the spatial coordinate ( leads to the line element 

1. The n > 4 case: Approach to homogeneous dust dynamics and black hole formation 

The n > A case corresponds to the black hole formation and the ratio FnP^ /\t\ approaches 
zero towards the singularity formation. From Eq. (13.161) . we have 

4 

^asym • (3.19) 

y 

In addition, from Eqs. (I3.8p - (I3.10I) . we obtain the solution for n > 4 as 

/o = ^ cosh^ (0 , i?o = ^ cosh^ (0 sinh (0 , (3.20) 

where / is an arbitrary constant and 

p = I sinh(C/3) X exp(-t/3), r^/^ = (Bf/A) cosh2(C/3) x exp(-2t/3). (3.21) 

Although the metric tensor written by Eqs. (I3.20p is derived as an asymptotic form in the 
limit t — > oo, it is identical with an exact solution of the Einstein equations describing 
spherically symmetric homogeneous dust collapse. This unfamiliar form of the metric tensor 
is due to the separable volume gauge; the timelike congruence parametrized by the coordi- 
nates {t,r) (or () covers only a limited region of the spacetime (see Fig. [1]). In fact, from 
Eq. fl3.14p . the radial velocity is = tanh((^/3) and this family of the timelike curves has 
the null boundary = 1 at C = oo (r ^ oo). 
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2. The n = 2 case: Growth of the central density gradient and naked singularity formation 



For n = 2, the ratio F2p'^/\t\ grows with the lapse of time and the inhomogeneity due 
to the term -F„p" is significantly involved in the central naked singularity formation. Under 
the approximation F2P^/\t\ ^ 1, we obtain the metric functions 

Iq = 1, Ro = lsmh{C) (3.22) 

and 

Fs/y/s = /sinhC X exp(-t), = (95/2/4) cosh^ C x exp(-2t). (3.23) 

This asymptotic form of Jq and Rq represents a fiat metric written in the accelerating 
coordinate system {t,(}- Using Eq. fl3.14p . we have = tanh^. The spacetime region 
covered by the accelerating coordinate system also has the null boundary at ( = oo. 

The timelike congruence parametrized by (t, () converges to the point r = p = with 
the four- velocity u = bq. The spatial distance between the two neighboring timelike curves 
( and C + '^C changes in proportion to exp{—t)d(. This change of the reference congruence 
may be larger than that of the matter contraction effect which increases the proper energy 
density fi. Thus, the effect of matter on the metric becomes weak and we have the fiat 
metric fl3.22p as the dominant asymptotic behavior. We can find from Eq. (13.1 6p that the 
Hubble normalized density Q asymptotically takes the form 

Q^^^ = ^ir-6/7;2/7^Qgj^2^ (sinhC)-^'/'e-2*/^ ^ . (3.24) 

This exponential decay of i^asym is a remarkable feature of the asymptotic evolution for the 
n = 2 case. As the contribution of matter to the metric is negligible, one can call the 
behavior represented by Eq. f l3.24p as the "vacuum-dominated" evolution |l6|. 

It must be noted, however, that the approximation F2p^ /\t\ ^ 1 breaks down near the 
center p = (C = 0) which is regular during r < (t < oo). Let us denote the asymptotic 
value of VL at the center as VLq. Then, we obtain VLq = 4/9 even for n = 2. Note that 
Eq. fl3.24p shows a infinite increase of Q in the limit ( 0. This increase of Q with respect 
to the spatial coordinate ( should be suppressed if the ratio F2p'^/\t\ becomes smaller than 
unity in the vicinity of the regular center. Unfortunately, it is difficult to see analytically the 
smooth decrease of Q from the central value Qq. Nevertheless, it is sure that the gradient of 
fl near the center increases infinitely as t increases. Hence, the exponential decay (I3.24p of 
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f2asym should be rather regarded as the growth of the density contrast between the central 
region near C = and the outer region ( ^ 1. We expect such a profile of fiasym with a large 
gradient with respect to the coordinate C is essential to the naked singularity formation. 



3. The n = 3 case: Existence of the critical density gradient between naked singularity forma- 
tion and black hole formation 

In this case, the ratio F^p^/lrl remains finite in the limit t oo. As r ~ ~ exp(— t) 
for this case, by introducing the ratio 

k = ^ (3.25) 

\t\ 

as a function of the spatial coordinate (, we obtain from Eq. (13.161) the asymptotic profile 
of Q as follows 

4 

^^^^"^^^^ = 9il + 3k)il + k) • ^^-^^^ 
By using Eq. ( 13. 14^ . the radial velocity v^' can be written as 

fcV3(i + 3A;) ^ F3 



26V3(i + A;)V3' 53/2- 
Then, using Eqs. (I3.8l) - (l3.10p . the equation determining k is given by 

dk k{l - {v''Y} 



(3.27) 



(3.28) 



Imposing the boundary condition = at C = 0, we obtain the solution k = k{() containing 
the parameter b. For n = 3 case, we also have \v^\ — > 1 in the limit C 00, and the value 
of k remains finite at C = 00. It is clear from Eq. (I3.26P that fiasym decreases monotonically 
as k increases. Let us denote the values of k and r2asym in the limit C ^ 00 as k^o and Qoo, 
respectively. These limiting values depend on the value of the parameter b. The profile of 
^asym as a fuuctiou of ( is shown in Fig. [2j 
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FIG. 2: The profile of riasym as a function of the spatial coordinate C. Each lines correspond to 
b < be (the solid line), b = be (the long-dotted line) and b > be (the short-dotted line). As the value 
of b increases, the density contrast 6 = Qq/Qoq increases. For 6 > 6c ^ 15, the resulting singularity 
becomes naked. 

Although Qasyra = 4/9 at the center C = irrespective of the parameter b, the value Q^o 
decreases as b increases. For 6 ^ 1, we have koo — 2\^/3 and fioo — 1/3&. Recall that b > be 
gives the condition for the naked singularity formation. The numerical calculation shows 
that this inequality corresponds to Qoo < — 0.03. Using the central value flo = 4/9 of 
^asym, the density contrast defined by the ratio 6 = Qo/Qqo has a critical value Sc- S < 6c 
corresponds to the black hole formation and 6 > 6c corresponds to the naked singularity 
formation. For 6 ^ 6c, the density contrast becomes very large and this leads to the naked 
singularity formation. This behavior is consistent with the asymptotic profile of Q discussed 
in the case n = 2.^ 

We finally comment on the self-similar evolution in dust collapse in terms of the Hubble 

^ Although the naked singularity formation occurs both for n = 2 and for n = 3 with the parameter b > be, 
its strength is known to be different in the sense of Tipler j2l|. Namely, the arising central singularity 
is called gravitationally weak for the former case but strong for the latter one Q . This difference of the 
strength can also be interpreted from the viewpoint of the asymptotic behavior of fi, which claims that 
the "vacuum-dominated" evolution under the condition 17 ^ around the central region results in the 
formation of "weak" naked singularity. 
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normalized variables. The spherically symmetric self-similar solutions are defined as solu- 
tions of the Einstein equations reduced to a set of the ordinary differential equations with 
respect to a dimensionless variable (see 22| for its details). For n = 3, the Hubble normal- 
ized density Q, as well as the other Hubble normalized variables {v^, A*", A^, E/, Eg^) is given 
as the function of the single coordinate r in the limit t ^ oo. Although we have derived the 
asymptotic behavior of and Q through the coordinate transformations fl3.8l) - fl3.10l) . it is 
possible to obtain all the asymptotic forms of these variables directly from an autonomous 
system of the evolution equations fl2.25l) - fl2.30p and the constraints fl2.31l) - fl2.33l) for 7 = 1. 
As S^'*" and q approach zero in the limit t 00, if these two variables in addition to variables 
with partial derivatives with respect to t are neglected in the equations, we arrive at a closed 
system of the ordinary differential equations for {Q,v^ , X"^ , , Eg^) using the dimensionless 
coordinate ({see the next section for the details). This closed system will be equivalent to 
the Einstein equations with the requirement of self-similarity.^ The result obtained in this 
subsection means that such a self-similar behavior should develop at the final stage t 00, 
even if the initial distribution function F{p) given by Eq. (13. 5p contains the higher order 
inhomogeneous terms with p"' for n' > 4 in addition to the term F^p^. 



D. Condition for the naked singularity formation 

Combining the result of our analysis for the asymptotic behavior of the Hubble normalized 
density function Q in the LTB solution, we propose the following conjecture for the condition 
of the naked singularity formation: 

The development of the large spatial gradient or density contrast in the asymptotic profile 
of the Hubble normalized density parameter Q in the separable volume gauge gives a sufficient 
condition of the naked singularity formation. 

We have confirmed this conjecture holds for the dust case. In the next section, we examine 
whether this conjecture can be extended to the gravitational collapse of perfect fiuid with 
pressure. 

^ In fact, we can check that the function r^asym given by Eq. (|3.26p has the same form as Eq. (3.13) in [3| 
except for the numerical factor. 
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IV. EFFECT OF PRESSURE ON THE HUBBLE NORMALIZED DENSITY 



In the previous section, we have discussed the asymptotic behavior of the Hubble normal- 
ized variables through the coordinate transformation from the LTB solution of dust collapse. 
Unfortunately, there does not exist such a generic analytical solution of perfect fluid collapse 
for the equation of state (12. 6p with 1 < 7 < 2. We therefore have to analyze directly the 
asymptotic solutions of an autonomous system of the evolution equations fl2.25l) - fl2.30l) and 
the constraints (l23TD - fl233D . 

As was schematically shown in Fig. [H under the separable volume gauge (12.241) for the line 
element (I2.2ip . we can assume the timelike congruence parametrized by the coordinates t 
and r converges to the central singularity arising at t = 00 and the possible asymptotic time 
dependence can be assumed to be J ~ J ~ i? ~ exp(— t). The Hubble normalized variables q 
and then go to zero in the limit t 00, while the other variables {E/, Eg^, A", X^, Q, v^) 
may become finite dependent only on r. We express their asymptotic forms as {E, A, A, fl, v) 
without the indices r and 6, and also introduce the coordinate ( given by Eq. (13.17^ . It is 
easy to obtain A = -3jnv/2G from Eq. fl^:^ and dE/dC = -{A + X)E from Eq. (^M>- 
The assumptions for the timelike congruence and the possible asymptotic time dependence 
of (J, J, R) allow us to arrive at the closed set of equations for {A, fl, v) as follows 

{(7-iy + l}K-(7-l)}^ = 

2^nv^A {27 - 3 - (7 - 1)^;'} - ^7(7 - 2)G-^n\ {l - (7 - 1^} (1 - v^) 

- 2nv {(37 - 4)(7 - 1) - (272 - 57 + 4:)v^} , (4.1) 

^{{^-iy + i}{v^-{^-i)} dv ^ 

l-v^ dC~ 
27(7 - l)GvA + ^7^(7 - 2)^^;^ - G {(7 - 1)(37 - 2) + (7 - 2^} , (4.2) 

^ = 3^G-^nvA - -G-i(7 - 1 + v^)n --n-A^ + 1 . (4.3) 
etc 2 2 

From Eq. (14. 2p . we can find that there exists the null boundary |f | = 1 of the timelike 
congruence in the limit —>■ 00 just in the same way as dust collapse. 

It should be noted that this set of equations for {A, Q, v) is derived under the assumption 
Q ^ 0. If we consider the case — in the limit t 00, we must use the equation 

E^ -A^ + 1 = (4.4) 
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instead of Eqs. (14.10 and (14.20 . The similar equation can be obtained for the case n = 2 
of the dust collapse. In fact, under the requirement f2 = A = 0, we have the solution 
E = —1/ sinh(^, A = coth(^ which corresponds to the fiat metric given by 03.221) . It was 
shown in the previous section that this "vacuum-dominated" evolution should break down 
in the vicinity of the regular center due to the large spatial gradient of fiasym- For the perfect 
fluid collapse, the effect of pressure should become important at least in the subsonic region 
f ^ < 7 — 1 near the center and works to suppress the development of the density contrast. 
Nevertheless, it is plausible that the "vacuum-dominated" evolution leading to the decay 

— is allowed in the supersonic region f ^ > 7 — 1, where pressure becomes ineffective. 
This means that a large gradient of fiasym generated in the supersonic region is in favor of 
the naked singularity formation. Unfortunately, it is very difficult to study in details the 
perfect fluid collapse allowing the decay of Vl with the lapse of time in the supersonic region 
and to check the validity of our expectation. Therefore, in the following, our discussion 
is restricted to the case that VL remains nonzero in the limit t— s>ooinO<C<oo, and 
Eqs. fliAD-fOD can be applied. 

The set of the ordinary differential equations f l4.ip -f H73l) with respect to the nondimen- 
sional coordinate C describes the self-similar behavior of the Hubble normalized variables 
such as Vt. Although in this paper, we assume that they become asymptotically dominant 
equations as a result of gravitational collapse started from general non self-similar initial 
conditions, they are equivalent to the exact self-similar Einstein equations which have been 
extensively studied in the comoving coordinate system. We can see the existence of a singu- 
lar point of the differential equations at f ^ = 7 — 1 where the coefficients of the derivatives 
dQ/d( and dv/d( in Eqs. (14.11) and (14. 2 p vanish. This corresponds to the sonic point where 
the fluid velocity relative to a tangent surface of a homothetic Killing vector equals to the 



sound speed [l7|, l23| , and the self-similar solutions may become singular. 



We consider the self-similar solutions regular in < C < 00 including the sonic point. 
The simplest example is the so-called flat Friedmann solution with 7 > 1. It is well known 
that the arising singularity becomes spacelike in the homogeneous perfect fluid collapse. For 
this solution, the asymptotic density function is given by 



4 

^asym = ^ 



l + 7sinh2{(37-2)C/37} 



(4.5) 



l + sinh2{(37-2)C/37} 
The ( dependence of this function is shown in Fig. 3. This form can be obtained through the 



20 



coordinate transformation from the comoving coordinate system to the coordinate system 
{t, r, 6, (f)} as was performed in Sec. IIIII Of course, it is also possible to derive this form 
directly from Eqs. fl4.ip - (14.31) with = — tanh{(37 — 2)(/3'~f}. This solution can be inter- 
preted as an extension of the case n > 4 of the dust collapse, for which we have fiasym = 4/9. 
By virtue of the effect of pressure, we have the density contrast 6 = Qq/Q^o = I/7 < 1 in 
this homogeneous perfect fluid collapse. 
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0.2 
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FIG. 3: The spatial profile of l^asym for the flat Friedmann solution with perfect fluid. Each Hnes 
correspond to 7 = 1 (the solid line), 7 = 3/2 (the dotted- line) and 7 = 2 (the short-dotted line). 

As the density contrast or the spatial gradient is not so large, the evolution cannot become 
"vacuum-dominated" and the resulting singularity is not naked. This example supports our 
conjecture. 

The effect of pressure becomes clearer if we consider the general relativistic Larson- 
Penston (GRLP) solution ITj. This is the numerically obtained self-similar solution which 
describes monotonic collapse of an inhomogeneous perfect fluid with the regularity imposed 
at the sonic point. This self-similar solution may correspond to the case = 3 of the dust 
collapse provided that 7 is treated as a free parameter instead of b. The central singularity of 
this solution becomes naked for 7 < 7c ~ 1.0105 (weak pressure), while it is hidden behind 
a horizon for 7 > 7c (strong pressure). 

To obtain i^asym corresponding to the GRLP solution, we numerically solve the ordinary 
differential equations (I4.ip - (I4.3I) by requiring the regularity of the solution both at the center 
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= and the sonic point f ^ = 7 — 1. The result is shown in Fig HI 



asym 
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FIG. 4: The profile of the Hubble normalized density Oasym corresponding to the GRLP solution. 
Each lines correspond to 7 > 7c (the solid line), 7 = 7^ (the long-dotted line), 7 < 7c (the 
short-dotted line). As the value of 7 decreases, the density contrast 5 = i7o/^oo increases. For 
5 > 24, the resulting singularity becomes naked. 



It is clear from Fig. H] that i^asym decreases monotonically as C increases, just in the same 
way as the case n = 3 of the dust model (see Fig. [2]). Further, we note that the density 
contrast 5 = Qq/Q^o in this self-similar perfect fluid collapse increases with the decrease of 
7. In particular, the critical value of 7 for the naked singularity formation is 7 = 7c — 1.0105 
and using this value, we can roughly estimate the corresponding critical density contrast to 
be 6c = Qq/Qoo — 24 ^ 1 {Qoo — 0.04). Recall that for the case n = 3 of dust collapse, we 
have the critical value 6c = Qq/Qoo — 15 ^ 1. Although the critical value of the density 
contrast will slightly depend on the property of collapsing matter, the asymptotic analysis 
of the perfect fluid collapse also supports our conjecture that the very steep change in the 
asymptotic spatial profile of Q characterizes the type of the arising singularity. 

The asymptotic behavior of the Hubble normalized variable with the separable volume 
gauge becomes self-similar provided that the asymptotic value of Q is not zero. This may 



be interesting in relation to the so-called self-similar hypothesis 



24( 1 which asserts that a 
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self-similar behavior should be dominant near the dense central region at the final stage of 
collapse starting from general initial conditions. 



V. SUMMARY AND DISCUSSION 

We have studied the asymptotic dynamics of the naked singularity formation in spher- 
ically symmetric gravitational collapse of perfect fiuid. We first have examined inhomoge- 
neous dust gravitational collapse described by the marginally bound LTB solution. As our 
main result, we have revealed the different asymptotic behavior of the Hubble normalized 
density fl depending on the type of the initial inhomogeneity. By comparing to the known 
causal structure of singularity arising in the dust collapse, we arrive at the important conjec- 
ture that the very steep decrease in the asymptotic spatial profile of Q is the characteristic 
of the naked singularity formation. The validity of this conjecture has been also supported 
in the perfect fiuid collapse with pressure. 

Let us remark on the "vacuum-dominated" case of the dust collapse (the case n = 2), for 
which Q goes to zero in the limit t — > oo and the arising singularity becomes naked. Because 
the central value Qq should remain nonzero, the vacuum-dominated evolution leading to the 
naked singularity formation has been interpreted as the development of the steep spatial 
gradient in the profile of Q with the lapse of time. For the dust collapse, as was mentioned 
in Sec. IIIIl the asymptotic behavior ^ with the naked singularity formation is rather 
generic. Unfortunately, for perfect fiuid collapse with pressure, the dynamical evolution 
from an initial state to the vacuum-dominated state — > is not confirmed. One may find 



that ri approaches zero only in the outer supersonic region \v^\ > y^^)^^^. The self-similar 
behavior may be then dominant in the inner subsonic region. To reveal the existence of 
the vacuum-dominated evolution and the structural change from the inner region to the 
outer one remains as an important problem to be investigated, especially, in relation to the 
criterion for the naked singularity formation. 

If no symmetry is assumed in gravitational collapse, one can assert that the generic 



singularity becomes spacelike from the so-called BKL conjecture [25|, |26| that states the 
Einstein equation becomes local near the singularity. In terms of the 1+3 orthonormal frame 
formalism considered here, the BKL conjecture has been reformulated in pjj and applied to 



various types of collapsing matter (e.g., the perfect fiuid with 7 = 2 |27|). Although the BKL 
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conjecture has been supported even in the vacuum-dominated case, the arising singularity 
is known to be null (see [3]) f*^^ spherically symmetric dust collapse if ^ 0. In this 
paper, we have claimed that the naked singularity is originated from the steep gradient in the 
profile of Q. For nonspherically symmetric collapse, of course, the freedom of gravitational 
waves becomes important in the vacuum-dominated case and the long-wavelength modes 
will be crucial to suppress such a steep density gradient and avoid the naked singularity 
formation. The BKL conjecture on the causal structure of singularity is based on the 
assumption that the spatial gradient of dynamical variables becomes negligible at the final 
stage of collapse. This may not be valid for some models of gravitational collapse including 
spherically symmetric models. 
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